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Abstract. We describe the formation of solitons in NJL-type models and dis- 
(N ■ cuss the influence of the regularization scheme on the stability of the solution. 

I We concentrate on models with non-local regulators in which stable solutions 

exist without introducing additional constraints. 

pi i' 1 Introduction 

D ■ 

i-^ ' The quark models of baryons can be divided into two main classes: (i) con- 

stituent quark models describing baryons as bound states of three massive 
quarks interacting via a phenomenological potential and (ii) chiral models with 
^ . , "bare" quarks surrounded by a cloud of either effective meson fields or qq pairs. 

' The chiral model distinguishes itself from the constituent quark models in that 

the baryon masses, as well as the constituent quark masses are generated dy- 
namically. The two classes of models describe hadrons at two different levels: 
at a higher level, the constituent quark model successfully predicts the spec- 
trum of baryons and their excited states; at a lower level, the chiral models are 
able to describe the vacuum, meson and baryon sectors with the same effective 
Lagrangian and they model the origin of the masses of constituent quarks and 
in particular the role the chiral mesons play in the constituent quark model. 

We present here the calculation of solitons in NJL-type models and dis- 
cuss how different regularization schemes influence the stability of solutions. 
We focus on the version with non-local quark interactions as suggested from 
the instanton-liquid model which supports stable solutions without introducing 
artificial constraints as in previous calculations ^. 
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2 Basic Properties of the NJL 

The NJL model assumes an attractive interaction between quarks constructed 
in terms of quark bilinears in a form that obeys the chiral symmetry. To describe 
low energy phenomena, it is enough to consider the non-strange quarks and 
keep only the scalar and pseudoscalar terms. (For a review of the model and 
its applications see ref. j^.) In the chiral limit, the Lagrangian becomes: 

Cnjl = qi^d^r) Q + J + iq^l5Taqf] , (1) 

where a = 1, ... 3 is the isospin index. The interaction is point-like and the UV 
divergences have to be removed by a suitable regularization prescription. 

To solve the model in the mean-field (Hartree) approximation it is practical 
to introduce auxiliary chiral fields 

S ^ -G qq and Pa = -G qijsTaq , (2) 

such that the Lagrangian takes the semi-bo sonized form: 

Cnjl = qP {idt -h)q-^ {S' + Pj) , (3) 

h=-ia-V + PiS + IJ^TaPa) ■ (4) 

The vacuum in the Hartree approximation is calculated by assuming a con- 
stant value, M, for (S) and evaluating (qq) in (H) by summing up the bubble 
graphs or, equivalently, by performing functional integration over the quark 
fields. This yields the so called "gap equation": 

^-24^^;^, (5) 

A 

where J2a denotes a regularized sum over quark states with a cut-off parameter 
A. Starting from a zero mass we have generated a finite quark mass. The 
model thus exhibits the spontaneous breaking of the chiral symmetry, yielding 
a nonzero value for the vacuum quark condensate (qq) — —M/G, and three 
Goldstone bosons (pions) corresponding to small oscillation of the pseudoscalar 
field around its vacuum value (Pa) = 0. The model predicts the pion decay 
constant, /^r, and for a finite current quark mass, m 7^ 0, also a finite pion 
mass satisfying the Gell-Mann-Oakes-Renner relation. Two of the three free 
parameters, G, A, or m can be eliminated by fitting and to^. The remaining 
free parameter is usually expressed in term of M, the so-called constituent 
mass. 

In calculations of mesons and baryons several regularization schemes have 
been used including a sharp cut-off in 3- or 4-dimensions, the proper-time and 
Pauli- Villars regularizations and non-local regulators as suggested from the 
instanton-liquid model. All of them require a rather low cut-off parameter, A ~ 
600 MeV, which indicates that it may be derivable from the underlying theory. 
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Indeed, the instanton-liquid model [Q predicts a non-local interaction between 
quarks with A related to the inter- instanton spacing p^^ ~ 600 MeV. The 
calculated vacuum properties (quark condensate, gluon condensate) suggest 
the value M should lie in the range between 300 and 400 MeV. 

The soliton, corresponding to a baryon, is constructed around three valence 
quarks that polarize the vacuum generating a spatially non-trivial configura- 
tion. The resulting mean-field potential may make the system stable against 
decaying into free quarks. This is not a mechanism of confinement but simply 
a mechanism to bind quarks. If this scenario is realized it means that the inter- 
actions of the quarks with the chiral field may play an equally important role 
in the formation of the baryon as the color confining forces. 

The existence of a stable solution is closely related to the choice of regular- 
ization scheme. A simple sharp cut-off does not yield stable solutions. Using the 
proper-time or Pauli-Villars regularization a stable soliton can be formed only 
by restricting the chiral fields to lie on the chiral circle, i.e. S{r)'^+Pa{r)'^ — M^, 
otherwise the scalar field acquires an arbitrarily low value at the origin pro- 
ducing a strongly localized state with a zero energy. Such a constraint does 
not have any justification in the model. It is only in the version with non-local 
regulators that solitons exist without introducing additional constraints. 

3 Solitons in the NJL Model with Non-Local Regulators 

In a non-local version of the NJL, the qq interaction is smeared by replacing 
the quark fields in the interaction part of ([^) by delocalized quark fields: 

^Pix) = {xm = {x\r\q) = ^ rfce"=-(fc|g) . (6) 

k 

Here is a regulator, diagonal in fc-space, and k is the Euclidean 4-momentum, 
k^ = u}^ + k'^. The instanton-liquid model predicts a form of the regulator. For 
model calculations it is simpler to take a Gaussian form, — e^^ 1"^^ . Using 
the Euclidean formulation the Dirac Hamiltonian takes the form 

-\a.-V ^rliSr + \-1^iiTarP<,r , r r (-9,2 - V^ yl) . (7) 

The key point is the construction of the valence orbit. Let us first study the 
free quark propagator. In the chiral limit (m = 0) the inverse propagator takes 
the form fcp7p + r{k'^)'^M. The pole (setting = 0, fcg = —M^) occurs at the 
solution of 

= M2g2M,V^^ . (8) 

The most striking feature of (^) is that the solution exists only below a critical 
value of M. Above this value, free on-shell quarks do not exist in the vacuum. 

The treatment of the quark propagator in a spatially non-homogeneous field 
configuration of the soliton is a non-trivial problem because of the presence of 
the time-dependent regulator in the Dirac Hamiltonian. For stationary chiral 
fields S and Pa the Dirac operator h is diagonal in the energy representation: 

h{u^) I ) = iL^ I ) . (9) 
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The lowest (valence) orbit is obtained as a solution of Eq. ^ for lu ^ Eq- 

^_Z+^e^(^°'+^')(5(r)+i75raP4r))eii^(^°'+^'' |Ao) = eo|Ao) . 

(10) 

The solution of Eq. (|l^) exists for all M beyond Mcr as shown in Fig. |^. 

The soliton is sought iteratively by assuming a hedgehog shape for the chiral 
fields. Starting from an initial guess for5(r) andF„(r) Eq. (|) is solved for the 
valence orbit, |Ao), and the sea orbits, \Xuj)- In the next iterations new values 
for the chiral fields are obtained using the Euler-Lagrange equations: 



G J 2n ^ \uj + e\ (w^) 



Pair) .J /X I I \ / I u \ , f duj ^ {>'u^\r\r)i/3j5Ta{r\r\X^) 
= NMMr\r).h.ra{r\r\X,) + j ■ 

Here zq is the residue factor: zq = (l — i deo(i^)/dti>|^^j^ J ^ . The convergence 
is obtained for all values of M beyond a critical value A/cr- 

In this approach the soliton acquires the correct baryon number Q . 



4 Properties of the Soliton 

Figure ^ (a) shows the soliton energy as a function of M for three regularization 
schemes: (i) the non-local model described in Sect. 3 yields stable solutions 
above M^r = 280 MeV which are also energetically stable since the energy of 
three free on-shell quarks is always higher than the soliton energy, (ii) the model 
with a non-local regulator of the Gaussian shape that depends only on the 3- 
momentum k possesses stable solutions above M k, 325 MeV; however, beyond 
M K, 600 MeV the solution is again unstable since it becomes energetically 
favorable for two quarks in the soliton to acquire high momenta, lower their 
effective mass, and escape from the third quark (Such an awkward behavior is 
a consequence of breaking the Lorentz invariancc and does not show up in the 
model with the 4- momentum regulator.) (iii) the model using the proper 
time regularization and the chiral-circle constraint has stable solutions above 
M = 340 MeV that are energetically stable only beyond M = 400 MeV. 

The model with the 4-momentum non-local regulator supports another type 
of solutions with one or two quarks polarizing the sea (Fig. |l|). The solution 
with three valence quarks is stable against disintegrating into solitons with a 
lower number of valence quarks. This does not hold for a soliton with four or 
more valence quarks since in such a case the additional quarks would have to 
fill the grand spin 1 orbit which is energetically less favorable. 

The calculation of observables in the model with non-local regulator is a 
non-trivial task since the Noether currents acquire additional terms originating 
from the momentum dependent regulator Fortunately, their contribution 
remains small for physically interesting values of M below 400 MeV. The calcu- 
lated properties reflect a relatively large size of the soliton compared to solutions 
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using e.g. the proper time regularization as well as to the experimental values. 
In our opinion this is a consequence of the mean-field treatment rather than a 
serious deficiency of the model. Any improvement of the approximation such 
as elimination of spurious center-of-mass motion or projection onto subspace 
with good spin and isospin may considerably reduce the size |q| . 



5 Implications for Other Models 



The solitons constructed in the chiral model with non-local regulators share 
several features found in more phenomenological models. From the model it is 
possible to derive a version of the linear cr-model that approximates well the 
full model [|. 

The effective quark mass (i.e the square root of + shown in Fig. ^) 
drops almost to in the center of the soliton for M between 300 MeV and 
350 MeV similarly as in models based on restoration of the chiral symmetry 
inside the baryon like different versions of the bag model. 

The Goldstone-boson exchange constituent quark model also finds a qual- 
itative support in the non-local model. The solitons consisting of only one 
valence quark (see Fig. ^ can be interpreted as constituent quarks. Above 
M « 310 MeV these solutions are dressed in a cloud of chiral mesons that 
generate attraction and a 3-quark soliton is formed. The interaction between 
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Figure 1. (a) The energy of the soliton for different regularization schemes plotted as 
functions of the parameter M using the 4-dimensional non-local Gaussian regulator 
(solid line), the 3-dimensional Gaussian regulator (dashed line), and the proper-time 
regularization (dashed-dotted line), (b) The energy per quark of the soliton with 
three valence quarks (solid line) and with one valence quark (dashed line) plotted as 
functions of the parameter M. Also shown are the corresponding valence energies. 
Mq is the free-space on-shell quark mass. All quantities in MeV. For each case the 
values of /tt and m,r have been fixed to their physical values. 



6 




Figure 2. (a) Self consistently determined fields and baryon densities (multiplied 
by 47rr^); (b) effective squared quark mass for M = 325 MeV (solid line) and M — 
750 MeV (dashed line), plotted as functions of the radial coordinate r. 



such objects calculated in the linear cr-model [g shows a typical behavior of 
the potential used in the Goldstone-meson exchange constituent quark model. 
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